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In the framework of the 1/N expansion we show that the Callan-Symanzik /3-function associated with the 
four-point coupling g is non-analytic at its zero, i.e. at the fixed-point value g* of g. This singular behavior can 
be interpreted by renormalization group arguments, and written in terms of scaling correction exponents. 

We obtain accurate determinations of g* in 3-d and 2-d by exploiting two alternative approaches: the e- 
expansion in the </> 4 formulation of the O(N) model, and the high-temperature expansion of the lattice iV-vector 
(O(N) nonlinear a) model. These results are compared with the available estimates by other approaches, such as 
the fixed-dimension perturbative expansion, Monte Carlo simulations, etc... 

We also present results for the n-point renormalized coupling constants that parameterize the behavior of the 
effective potential in the high- and low-temperature phases. 
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The renormalization-group theory of critical 
phenomena provides a description of statistical 
models in the neighbourhood of the critical point. 
For 0(N) models calculations are based on the 
(/> 4 -field theory. A strategy, which has been 
largely employed in the study of the symmetric 
phase, relies on a perturbative expansion in pow- 
ers of the zero-momentum four-point coupling 
g performed at fixed dimension. The theory is 
renormalized at zero momentum by requiring 

T (2) (p) a0 = 8 a0 Ze 1 [m 2 + P 2 + 0(p i )] (1) 
r< 4 > (0,0,0, 0) aPlS = Z G 2 mg S a0lS /3 (2) 

where Sap^s — Sapd^s + S a7 Sps + 5 a g5py. When 
in — ► the coupling g is driven toward an infrared 
stable zero g* of the Callan-Symanzik /3-function 

/3(g) = mdg/dm\ g0t A- (3) 
g* is also obtained as the critical limit of 
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where x 1S the magnetic susceptibility, 
second-moment correlation length, and xa = 
Si,j,k(^(°) ' (t>{xi)(l>{xj) ■ </>(%k))c- We recall that 
N- vector (nonlinear O(N) a) and <fi A models de- 
scribe the same critical behavior. 

In the framework of the 1/N expansion the 
analysis of the next-to-leading order shows that 



the Callan-Symanzik /3-function is non-analytic 
at its zero, i.e. at the fixed-point value g* of 
g 0] . The large- N result agrees with the singular 
behavior 



m = 

+ci{g* 



-v{9* ~ 9) 
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analytic terms 
+ d 1 (g*-g)^ +.(5) 



(A = wv and A2 are scaling correction expo- 
nents) that can be derived using renormalization 
group arguments 

A precise determination of g* is crucial in 
the field-theoretic approach based on the g- 
expansion, where the critical exponents are ob- 
tained by evaluating appropriate (resummed) 
anomalous dimensions at g* . In this approach 
the resummation of the g-expansion is usually 
performed following the Le Guillou Zinn-Justin 
(LZ) procedure |3), which assumes the analytic- 
ity of the /3-function. The presence of confluent 
singularities may then cause a slow convergence 
to the correct fixed-point value, leading to an un- 
derestimation of the uncertainty. A more gen- 
eral analysis explicitly allowing for the presence 
of confluent singularities would slightly change 
the value of g* for small values of N and conse- 
quently the values of the critical exponents . It 
is therefore important to exploit other approaches 
to the study of O(N) models, which can provide 
a check of the estimates of g* from the resumma- 
tions of the ^-expansion. We considered two al- 
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Table 1 



Three-dimensional estimates of g* = g*(N + 8)/(487r). For the Ising model we also mention the recent 
Monte Carlo estimate g* — 1.397(14) (actually derived by us using the data kindly made available by 
the authors). A more complete list of the available estimates of g* can be found in Ref. (IJ. 
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3-exp.|4 
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H.T. cubic|5] 


H.T. bcc|5] 
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2 
3 
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1.390(17) 

1.397(8) 

1.413(13) 

1.387(7) 

1.366(15) 


1.413(6) 
1.411(4) 
1.403(3) 
1.391(4) 
1.377(5) 


1.39 
1.40 
1.40 
1.39 


1.393(20) 

1.406(9) 

1.415(11) 

1.411(12) 

1.396(16) 


1.388(5) 

1.408(7) 

1.411(8) 

1.409(10) 

1.392(10) 


1.387(5) 
1.407(6) 
1.406(8) 
1.406(8) 
1.394(10) 


Table 2 

Two-dimensional estimates of g* = g*(N 
asymptotically free, thus j3 c = oo. 


+ 8)/(24vr). We 


recall that 0(N) a models with N > 3 are 
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e-exp. 




g-exp. 
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M.C. 


1/N-exp. 
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3 
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1.69(7) 
1.75(5) 
1.79(3) 
1.72(2) 
1.64(2) 


1.85(10) | 
1.749(16) D 


1.679(3) | 
1.7540(2) [ 
1.810(10)1 
1.724(9) | 
1.655(16) [ 
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1.71(12) ra 

1.76(3) p| 
1.73(3) |l| 


1.758 
1.698 



ternative approaches: the e-expansion in the con- 
tinuum <^> 4 formulation and the high-temperature 
(HT) expansion of the AT- vector (lattice 0(N) a) 
model. We extended the e-expansion of g* to 
0(e 4 ). Accurate estimates of g* in 3-d and 2-d 
were obtained by a constrained analysis of the e- 
series using its known values at lower dimensions. 
Moreover we reanalyzed the available HT expan- 
sion of g a in the Af-vector models, by a method 
able to handle the leading confluent singularity 
(for a more recent analysis using longer series see 
Ref. |J) . In Table [j] and || we present our 3-d and 
2-d results respectively. For comparison we also 
report some of the available estimates from other 
approaches. The agreement among the various 
estimates of g* is globally good. 

The results in Table [l] indicate that the system- 
atic error in the LZ resummation due to the non- 
analytic terms in Eq. (||) should be small. This 
may be explained by the fact that, for small val- 
ues of N, the exponents in Eq. (||) are close to in- 
teger numbers, indeed A 2 /A~2, A 3 /A~3, and 
1 + 1/A ~ 3. However, the results for N = 0, 1 
are slightly lower than the estimates given by the 
LZ resummation of the 0(g 7 ) series of (3(g), thus 
favouring the more general analysis of Ref. H]. 



This would lead to a small change in the esti- 
mates of the critical exponents. For instance, in 
the case N = (self avoiding walks) the resum- 
mation of the 0(g 7 ) series of 7(g) evaluated at 
g* = 1.413(6) gives 7(5*) ~ 1.160. A lower value 
g* ~ 1.39, as indicated by our calculations, would 
lead to 7(3* = 1.39) ~ 1.158, in substantial agree- 
ment with the recent result of Monte Carlo sim- 
ulations 7 = 1.1575(6) 12 and with the analysis 
of the e-expansion: 7 ~ 1.158. 

The effective potential is widely used in the 
field-theoretic description of fundamental interac- 
tions and phase transitions. In statistical physics 
it represents the free-energy density J 7 as a func- 
tion of the order parameter. T can be expanded 
in powers of the renormalized magnetization tp 
1 
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where AJ 7 = !F((f) — !F(0), and g%j are the zero- 
momentum 2j-point renormalized coupling. By 
definition g% = \ and gi — g- Setting ip = 



,(d-2)/2. 



/ y/g and r2j = g2j/g J we write 
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Table 3 



Three-dimensional estimates of r@ and rg. A more complete list of the estimates of r-ij is reported in 
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e-exp. 


»'6 

g-exp. 


ERG 


H.T. 


e-exp. 


rs 
g-exp. 


ERG 


H.T. 
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2.058(11) 


2.053(8) 


2.064(36) 


1.99(6) 


2.48(28) 


2.47(25) 


2.47(5) 


2.7(4) 


2 


1.94(11) 


1.967 


1.83 


2.2(6) 


3.5(1.3) 




1.4 




3 


1.84(9) 


1.880 


1.74 


2.1(6) 


2.1(1.0) 




0.84 




4 


1.75(7) 


1.803 


1.65 


1.9(6) 


1.2(1.0) 




0.33 





In order to evaluate the first few r2j we per- 
formed a constrained analysis of the e-expansion 
of r2j ■ In Table || we report our 3-d results for 
tq and rs- We compare them with some of the 
available estimates from other approaches, such 
as d = 3 g-expansion [p|fl4|], approximate solu- 
tion of the exact renormalization group equation 
(ERG) (l^JT^], high-temperature expansion [|l7[ 
|l8| . In the case of the Ising model also rio has 
been roughly estimated: t\q — —20(15). 

In two dimensions an analysis of the HT expan- 
sion of the free-energy of the Ising model in the 
presence of an external field gave r 6 = 3.678(2), 
r 8 = 26.0(2) and r w = 275(15). These num- 
bers compare well with the estimates r 6 = 3.69(4) 
and r§ = 26.4(1.0) obtained from our constrained 
analysis of the e-expansion. Moreover we ob- 
tained r 6 = 3.54(7) and r 8 = 25.1(2.0) for 
N = 2 (2-d XY model), and r 6 = 3.33(6) and 
r 8 = 20.3(1.7) for N = 3. 

In the broken phase of the 3-d Ising model the 
effective potential at the coexistence curve can be 
expanded as 



i=2 



J 



g-(v-V y(.8) 



where <7 ■ are the zero-momentum j-point renor- 
malized coupling in the broken phase (g^ = 1 
by definition) . A constrained analysis of the e- 
expansion 19 gave the estimates g$ = 13.06(12) 
and g~l = 75(7). The parametrization (||) does 
not apply to the case N ^ 1, due to the presence 
of Goldstone bosons. One indeed finds 

d/(d-2) 



?{<P) ~ -^(Vo) ~ c (ip 2 - tpl) 



(9) 



In 3-d, where d/ (d—2) = 3, logarithms appear in 
the corrections to the leading behavior. 
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